Odd-even staggering of binding energies as a consequence of pairing and mean-field 

effects 
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Odd-even staggering of binding energies is studied in finite fermion systems with pairing correlations. 
We discuss contributions of the pairing and mean-field to the staggering, and we construct the 
binding-energy filters which measure the magnitude of pairing correlations and the effective single- 
particle spacings in a given system The analysis is based on studying several exactly-solvable many- 
body Hamiltonians as well as on the analytical formulas that can be applied in the weak and strong 
pairing limits. 

PACS number(s): 21.10.Dr, 21.10.Pc, 21.60. Jz, 71.15.Mb 



I. INTRODUCTION 

The odd-even staggering (OES) of binding energies, 
which reflects stronger binding of even-particle-number 
systems than their odd-particle-number neighbors, has 
long been known in atomic nuclei W and recently has 
been observed also in metal clusters 1^ and in ultrasmall 
superconducting grains However, as discussed re- 
cently the apparent similarity between the OES in 
all these finite many-fermion systems is deceptive. Al- 
though two basic physical mechanisms are involved: {i) 
an effect of spontaneous breaking of the spherical sym- 
metry (Jahn- Teller effect |^) and, [ii) blocking of pair 
correlations by an unpaired fermion, the origin of the 
OES is system-specific. 

In small metal clusters, OES is believed to be mainly 
due to the non-spherical shape of the underlying mean 
field 1^-^ and, thus far, neither the empirical evidence 
nor the theoretical calculations j|] support the presence 
of superconductive correlations in metal clusters. In su- 
perconducting grains, on the other hand, OES is believed 
to be predominantly due to a blocking effect caused by 
the presence of an odd electron (see Refs. lTo|-[l7||). In 
the smallest grains, where the average spacing d of the 
electronic energy levels becomes comparable to the size of 
the pairing gap [|l8|, d ^ A, the OES must invoke both 
mean- field and pairing effects. In the case of metallic 
grains, it is, however, difficult (if possible at all) to pin 
down the detailed structure of the single-electron level 
spectrum. Most of the theoretical models applied to the 
problem of superconductivity in grains use an equidis- 
tant level spectrum because of its analytical simplicity 
(see Refs. |l^Jll|,|l^,^ and references quoted therein). 
However, more realistic treatments require a non-uniform 
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distribution of single-particle levels as is done, for exam- 
ple, in Ref. |l| (see also Ref. |§). 

A slightly different situation holds in atomic nuclei 
due to the presence of two types of fermions. The 
strong and attractive effective proton-neutron interac- 
tion gives rise to the appreciable symmetry energy con- 
tribution [^{N —Z^] to OES. Because the symmetry en- 
ergy smoothly varies with particle number and is almost 
insensitive to shell effects |^, the standard way of ex- 
tracting the OES has been by means of the higher-order 
binding-energy indicators (filters), such as the four-point 
expression of Refs. ||2l|,|2^ (see also Ref. Q and ref- 
erences therein). This classical reasoning, which stems 
from liquid-drop or Fermi-gas models, has recently been 
questioned Namely, it has been demonstrated, using 
fully self-consistent Hartree-Fock theory, that the con- 
tribution to OES due to symmetry energy is, in fact, 
nearly cancelled out by the contribution coming from the 
smoothed single-particle energy. 

The aim of this work is to demonstrate the consistency 
of our new interpretation using the results of simple, ex- 
actly solvable models. Section || discusses the binding- 
energy indicators used in this study. It also presents 
the systematics of experimental pairing gaps obtained 
using these filters. The many-body models employed in 
this work (seniority model, equidistant-level model, and 
pairing-plus-quadrupole model) are described in Sec. Ill , 
together with the results of the calculations. A short 
summary and conclusions are given in Sec. 0. 



II. ODD-EVEN STAGGERING FILTERS 

The simplest way to quantify the OES of binding en- 
ergies is to use the following three-point indicator: 

A(3)(Ar) EE "^[BiN - 1) + B{N + 1) - 2B{N)l (1) 
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where ttn = (—1)^ is the number parity and B{N) is the 
(negative) binding energy of a system with N particles. 
In expression (|^), the number of protons Z is fixed, and N 
denotes the number of neutrons; i.e., this indicator gives 
the neutron OES. An analogous proton OES indicator is 
obtained by fixing the neutron number N and replacing 
iVbyZin®. 

By applying filter (|l|) to experimental nuclear bind- 
ing energies, -Boxp , one obtains the experimental neutron 

and proton OES, Acxp. Similarly, by applying this filter 
to calculated binding energies, one obtains information 
about theoretical results for the OES. In the following, 
the same subscripts are used to denote the binding en- 
ergies and values of indicators (j^) obtained for a given 
model. For example, for the single-particle model de- 

(3) 

scribed below, the resulting values are Bsp and Asp . 
Note that by using filter (Q) we only aim at facilitat- 
ing the comparison of calculations with data; however, 
in essence we always compare and analyze experimental 
and calculated binding energies. 

As a simple exercise, let us first calculate the OES for 
a system of particles moving independently in a fixed de- 
formed potential well. In such an extreme single-particle 
model, the binding energy is 



k^k, 



(2) 



fe=i 



where (=0, 0.5, or I) and stand, respectively, for 
the occupation number and single-particle energy of the 
/c-th, two-fold degenerate level, and the particle number 
N is given by 



N : 



(3) 



fe=i 



In the above expressions, the single-particle en- 
ergies Ck {k—1,2,...,^) appear in ascending order 
(ei<e2,. . . ,<esi), and denotes the number of twofold 
degenerate single-particle levels. 

For the ground state of an even-A^ system the N /2 
lowest levels are filled. In the neighboring odd- A*" system, 
the odd particle occupies the lowest available level. This 
implies: 



Ag)(Ar = 2n + f) = 
Ag)(7V = 2n)-i(e, 



(4a) 
(4b) 



Hence, in the absence of the two-body interaction, filter 
(|l|) vanishes at odd particle numbers, and it gives half 
of the single-particle level spacings at even particle num- 
bers. Therefore, following Ref. [Q, we assume that any 
other effect leading to the OES beyond the pure single- 
particle (mean-field) contribution is characterized by val- 
ues of (|^) at odd particle numbers. We also assume that 
such an effect varies smoothly with the particle number, 
and, therefore, by subtracting filter (M at odd and even 



particle numbers, one may obtain information about the 
single-particle level spacing at the Fermi level. This hy- 
pothesis will be verified in Sec. [II using several theoret- 
ical models. 

Consequently, at odd particle numbers iV = 2n-|- 1, we 
shall use the following filter, 



A{N) = A^^Hn), 
while at even particle numbers N = 2n we define. 



(5) 



= 2A(3)(iV)-2A(3)(iV±l 



(6a) 
(6b) 



6e{N) = 2A(3)(7V) - A(3)(7V - 1) - A^^\N + 1). 

Filters Se^ correspond to subtracting values of filter (|^) 
at particle numbers higher or lower by one, while Se em- 
ploys the symmetric average of both. The differences 
between these three filters reflect, therefore, small vari- 
ations of the OES effects with particle numbers. Up to 
these small variations, we associate the values of these 
filters with the spacings of the single-particle levels. 



de^{N) ~ de{N) ^ ejv/2+1 - 6^/2- 



(7) 



It is clear that filters Se^ use masses of four nuclides 
near the given particle number N; i.e., they constitute 
asymmetric expressions in which either heavier or lighter 
nuclides dominate. On the other hand, Se uses masses 
of five nuclides symmetrically on both sides of N. The 
advantages of using either of these filters depend, there- 
fore, predominantly on the availability of experimental 
data in the isotopic or isotonic chains. For instance, in 
Ref. dealing with light- and medium-mass nuclei, the 
filter 6e~^ was discussed. 

It is instructive to relate the asymmetric energy- 
spacing filters 5e^ of Eq. ( |6b| ) to the particle separation 
energies SiN)=B{N-l)-B{N), i.e.. 



6e+{N) = S{N) - S{N + 2), 
6e-{N) ^ S{N - I) - S{N + I). 



(8a) 
(8b) 



One sees that the filter 6e^ depends on separation ener- 
gies of particles from even systems, while the filter Se~ 
depends on those from odd systems. This does not in- 
volve any asymmetry in treating even and odd systems, 
because, obviously, every particle-separation energy de- 
pends on one mass of an even system and on one mass 
of an odd system. 

By using the three-mass filter (||), we hope to avoid 
mixing the contributions to the OES which originate from 
single-particle structure with those having other roots, 
e.g., pairing correlations. Moreover, because it involves 
three masses only, this filter allows obtaining experimen- 
tal information on the longest isotopic or isotonic chains. 
Figures |l] and || display experimental values of the neu- 
tron A^^p and proton A^^p OES (|). (It is to be noted 
that these results slightly differ from those presented in 
Ref. M . Firstly, experimental masses were taken from an 
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updated mass evaluation . Secondly, only the masses 
having an experimental uncertainty less than 100 keV 
have been considered.) 

It is seen that, especially for the light- and medium- 
mass nuclei, there exists a substantial spread of results 
around the average trend. This suggests that neutron 
and proton OES effects are not only functions of neutron 
and proton numbers, respectively, but that a significant 
cross-talk between both types of nucleons exists. Numer- 
ous studies of this isotopic dependence of the OES exist 
p^-p6[, usually based on higher-order filters such as the 
four-point mass formula pl|j2^ , |2^ : 

A(4)(iV) EE ^ [WiN - 1) - 3B{N) - B{N - 2) 

+ B{N + 1)] = i[A(3)(iV) + A(3)(iV - 1)]. (9) 

Unfortunately, as demonstrated in Ref. Q, the higher- 
order filters mix the pairing and single-particle contri- 
butions to the OES; hence, they are not very useful for 
the purpose of extracting the pairing component. Since 
the detailed analysis of the isotopic dependence of the 
three-mass filter results is not yet available, below we 
discuss the average trends; i.e., for each N the values of 
AJ^xp(A^) are averaged over all isotones for which the data 
exist, and, similarly, for each Z the values of AJ^p(Z) are 
averaged over all isotopes. 

Figure ^ (and solid lines in Figs. |] and ||) present such 
average values, Acxp, for the neutron and proton OES. 
Neutron and proton values of the OES follow a similar 
pattern. Namely, they systematically decrease with A, 
and they are reduced around shell and subshell closures, 
as expected. It is also seen that neutron pairing gaps are 
systematically larger than the proton ones. The shift is 
mainly due to a mass dependence of the pairing strength 
and results from the fact that at a given value N , data 
are more available for lighter nuclei than at an identi- 
cal value of Z. Indeed, for N^Z, data exist mostly at 
A{N)<A{Z). A weak contribution from the Coulomb 
energy is also expected to contribute to this shift. 



In Eq. (11) G is the pairing strength parameter. 



(12) 



fc=i 



denotes the monopole-pair creation operator, and k de- 
notes the time-reversed state. 

Properties of the Hamiltonian (nm depend on the ratio 



G 



(13) 



where k represents the strength of Hq. For both 77 ^ 
1 (weak pairing) and 77 ^ 1 (strong pairing), one can 
treat the Hamiltonian ( p^ perturbatively. However, the 
situation encountered most often in the nuclear physics 
context is the intermediate case (7;'^0.4) in which pairing 
correlations are strongly infiuenced by the nuclear mean 
field. 



A. Limiting cases 

To facilitate the discussion of results of exactly solvable 
models, it is instructive to consider weak- and strong- 
pairing limits of (0). To this end, we introduce the ex- 
plicit single-particle Hamiltonian 



(14) 



fc=i 



where are the single-particle energies defined in Sec. ||. 



1. Weak-pairing limit, 77 >C 1 



In this limit, the obvious expansion parameter is rj. 
The unperturbed ground-state energy, B^^^^ of an even- 
even system with particle number N—2n [n stands for 
the number of pairs) is that of Eq. (^), i.e., 



III. THEORETICAL MODELS 

In this section, we investigate several exactly solvable 
models to see the interplay between the particle-hole and 
particle-particle channels of interaction. In all cases, the 
Hamiltonian has the form 



H — Hq + Hr, 



(10) 



(15) 



The pairing Hamiltonian scatters the nucleonic pairs 
from hole states (7,7) to particle states (j,j), and the 
corresponding scattering matrix element is constant and 
equal to —G. The first order of perturbation theory gives 
the energy correction 



where H = Hq is either the intrinsic (i.e., deformed) 
single-particle Hamiltonian or the laboratory-system 
quadrupole-quadrupole Hamiltonian, and -ffpair is always 
the monopole-pairing (seniority) Hamiltonian: 



H ■ = -GP^ P 



(11) 



while the second-order correction is 



B 



(2) _ 



E 



(16) 



(17) 



j—n-\~l i—1 
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In the case of an odd system with N=2n+1, the anal- 
ogous expressions can be obtained by the following sim- 
ple modifications. First, in the zero order, the (n+1)- 
th level is occupied (blocked) by one particle, and hence 
the single-particle energy e„+i should be added to B^^^. 
Moreover, since the pairing Hamiltonian docs not couple 
orbitals occupied by one nucleon (the blocking effect), 
the orbital containing the odd particle must be excluded 
from the sum in Eq. (p^, and the number of pairs in 
Eqs. (|l|) and (0) must become n=(iV-l)/2. 

Adding together the zero-, first-, and second-order con- 
tributions to the binding energy, one obtains the corre- 
sponding expressions for A^'^h 



1 

J— n+2 



E 



i=l 



(18a) 



A(f^),(2n) = ^!i±Y^ + iG 



E^ — + E — 



j=ri+l 



(18b) 



Finally, for the single-particle energy-splitting filters 
( |6a| ) and (6b), one obtains: 



6e 



>)<i 



- e„+i 



2 



n ^ n 
E g _ e- E 



j=n+l 



G2 



+- 



1 

E — ^ 



n— 1 



G2 
" 4 



E^^- E 



fe=l 



efc - e 



Cfc - e„+i 



(19a) 



(19b) 



(19c) 



It is clear that in the limit of small G (small compared 
to the typical single-particle energy spacing) , the energy- 
spacing filters correctly extract the single-particle spec- 
trum from the total binding energies. 



Strong-pairing limit, 3> 1 

In this case, the expansion parameter is 77^^. In the 
zero order, the wave function of an even system with 
N—2n is the ground state of Hpair', i-e., it corresponds 
to the state with the maximal quasispin L,„ax=^^/2 and 
the third component of quasispin LQ={2n~D,)/2. The 



binding energy is given by the se niority-model expression 
(see Refs. [|^|| and Sec. |IIIB| ), 



4>^i(A^ = 2n) = -G[L„,ax(i 



max V-^ max 



+ l)-io(io-l)]. 

(20) 



In the quasispin formalism, the single-particle Hamil- 
tonian is a combination of a scalar and a vector operator 
(with respect to the quasispin group), and this implies 
the AL=0,±1 selection rule for its matrix elements. The 
expectation value of Hsp in the lowest L=Linax state is: 



B. 



(1) 



where 



2n) = (Lmaxio|^sp|imaxio) = 2ne, (21) 



(22) 



k=l 



is the average single-particle energy. The second-order 
correction to the energy is given by 



Bl^l.iN = 2n) = - 



E 



I {LvaaxLo 



L<L„ 



E 



LLo 



(23) 



In Eq. ( P2|) a denotes all the remaining quantum num- 
bers other than L and Lq, and E'^j^^ is the unperturbed 

excitation energy of \LLoa). Since Hsp can only connect 
the seniority-zero ground state with the seniority-two, 
i=imax-l states at energy E* = Gil, Eq. ( p3| ) can be 
reduced to a simple form; 

b1^1,{N = 2n) - -^(L„,axio|i?sp|imaxio) 



Gn 
1 

Gn 



(-^maxiol^^spl-^maxio)^- (24) 



By introducing the variance of single-particle levels. 



n 

k=l 



(cfc 



(25) 



one can derive a simple expression for the second-order 
correction: 



B. 



(2) 
J)»l 



{N^2n) = -^4n{n-n)al (26) 



That is, the first- and second-order corrections to the 
binding energy are given by the first and second mo- 
ments of the single-particle energy distribution. Note 
that for the degenerate j-shell, the second-order correc- 
tion is zero, as expected. 

For odd particle numbers, A^=2n-f 1, one should take 
Lmax=i^-l)/2 and Lo={2n + l-n)/2 in Eq. (|^). As- 
suming that the odd particle occupies level n+1, this 
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level is removed from the sum of Eqs. ( |22D and (^5|), i.e., 
one needs to consider the remaining — 1 levels only: 



n - 



n 

- S 



1 



efe 



fc=i 

k^n+l 



n - 1 



(e-e„+i) (27) 



and 



k=l 
1 



n-i 



(e - e„+i)^ 



(28) 



The resulting corrections to the binding energies can be 
written as 



Bl^l^iN = 2n + 1) = 2ne' + e„+i 



(29) 



and 



G(f7-1)2 



4n{^ -n - ly- 



(30) 



The zero-order expressions for A*^"^) and Se^ in the 
strong pairing limit are given in Sec. Ill B| . By adding 
the zero- and first-order contributions to the binding en- 
ergy, one obtains the strong-pairing-limit expressions for 

A(3). 

A^^l, (2n + 1) = 1g» + " ~ ^ (e„+i - e) (31a) 



Q — 2n ( e. 



n - 1 



-e| :77^^ 7^. (31b) 



2{n - 1) 



It can easily be shown that the first-order correction to Se 
vanishes. Consequently, the seniority-model expressions 
discussed below give a good approximation to the single- 
particle energy-splitting filters ( |6a|) and (6b). 



B. Degenerate shell: Seniority model 

Let us consider the seniority (or pairing quasispin) 
model ll^,^, i.e., the model for N nucleons moving 
in a 2il-fold degenerate shell described by the seniority- 
pairing Hamiltonian (11). For this model the exact solu- 
tion can be written in terms of quasispin quantum num- 
bers [see Eq. (po|)]. Alternatively, the ground-state en- 
ergy can be expressed in terms of the seniority quantum 
number s (see, e.g., psf p. 222): 



B,{N) - ^\G{N - s)i2n -s-N + 2), 
s = for N ^ 2n, 



where 



s = l for N = 2n+l. 



(32) 



The corresponding value of A'^'^^ is given by: 



Af)(7V) = 



^GQ + iG for N=2n, 



iGn 



for N=2n+1, 



which implies 



6ef{N)^6e,{N)^G. 



(33) 



(34a) 
(34b) 



It is seen that, in t he se niority model, filters ^ and (||) 
give the OES values ( |34a| ) and single-particle energy spac- 
ings ( |34b| ) which are in depe ndent of the particle number 
N. Values of the OES ( |34a|) correctly reproduce the ex- 
act pairing gap A'=^Gn, which is defined as a half of the 
lowest excitation energy in an even system, and which we 
denote by a prime to distinguish it from the OES. The 
meaning of Scs is less obvious. It is because the mean- 
field (Hartree-Fock, HF) treatment of Hamiltonian ( pT| ) 
yields only one 2ri-fold degenerate single-particle level 
at energy — G, while our interpretation of Scg assumes 
that only the Kramers degeneracy is present. Neverthe- 
less, one may compare exact values of Scs with those 
obtained in an approximated way and see whether the 
approximate ground-state energies reproduce features of 
the mass spectrum represented by Ag and Seg- 

In the seniority model, the Hartree-Fock-Bogoliubov 
(HFB) equations (which in this case are identical with 
the BCS equations) can be solved analytically. Indeed, 
the BCS occupation coefficient is given by 



N -s 



2{n - s) 

(see e.g. Q p. 233), and the ground-state energy is: 



(35) 



s-l-BCS 



(N) 



(36) 



The resulting three-mass filter (I) can be written as 



hGn 



kGn 



for iV=2ri, 



f for N=2n+\, 



and hence filters (||) and (^) give 

{n - 1)2 + Nf 



A,+Bcs(A^) = ^Gf! 

1 

(5e,+Bcs(iV) = g\i 



n{n - 1) 



(37) 

(38a) 
(38b) 
(38c) 
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It is seen that i n th e Hmit of large (fi^l), the BCS 
approximation (38a) rep rodu ces the leading order 0{^l^) 
of the exact OES result (|34|). The 0{fl^) deviation from 
the exact result smoothly depends on N and reaches the 
minimum (~G/4) at the middle of the shell (N—fl). 

For the energy-spacing filters the deviations behave as 
©(ri^^). While the symmetric filter ^Cs+bcs does not 
depend on iV, both (Se^g^g and Je^g^g vary weakly 
with the particle number. Namely, at the beginning 
of the sheU (iVwO), ^e^+gcs^G and Se^+scs-G - 
at the middle of the shell def_^j^^.g—Ses+Bcs^ and at the 
top of the sheU (A^«217) 6e+_^^^^KiG and 6e^^^,,^KiG - 
2G/Vt. This behavior follows from a simple identity 
<5e^j3cs(^)— '^6g^_3f-,g(217 — N) reflecting the particle- hole 
symmetry of the model. 

The analysis presented in this section illustrates the 
advantages of comparing exact and approximate results 
(similarly as experimental and theoretical results) by 
looking at appropriate filters. Analytical results avail- 
able in this model allow the explicit study of pairing and 
mean-field effects. The high degeneracy of the seniority 
model does not allow for extracting the energy spacings 
between the deformed levels; to this end, results for the 
non-degenerate models are shown in the following sec- 
tions. 



C. Equidistant- level model (infinitely many levels) 

Let us consider a phase space of infinitely many, 
doubly-degenerate, equidistant single-particle levels 
spaced by d. Suppose that all the levels up to a certain 
Fermi energy are occupied by the fermions^ and that they 
interact through a two-body interaction V . The Hamil- 
tonian of the model reads. 



H ^2 



E 

A:— — oo 



(39) 



where N-^ is the number operator of the n-th level. Fur- 
thermore we assume that V is invariant with respect to 
the shift in single-particle indices: 



— ^fei+fc,fc2-l-fc;fc'i+fc,fe2+*:' 



(40) 



where fc is an integer number. For a moment we do not 
need to specify this interaction; let us only remark that 
the standard seniority-pairing interaction ( [ll| ) obviously 
obeys the shift-symmetry condition. 

Of course, we do not intend to consider here infinite 
numbers of particles and infinite energies, cf. Eq. (|39|). In 
practice, we should assume that the number of fermions 
interacting through V is finite but large. The remaining 
ones (e.g., occupying the most bound shells) form an inert 
core. This guarantees the effect of the finite spectrum not 
to be important. Consequently, in future discussions we 
assume that A^^l and that the first level belonging to 
the "interaction-active" space has k—\. 



Irrespective of how complicated interaction V is, its 
assumed shift-symmetry allows for an exact analysis in 
terms of filters based on ground-state energy differences. 
Indeed, for any particle number TV, the ground-state en- 
ergy reads 



Bieq(7V) = 2 kdNk{N) + V{N), 



(41) 



fc=-c 



where iV"fc(7V) and V{N) are the average values of particle 
number operators Nk and interaction V in the ground- 
state wave functions for N particles. Since the particles 
interact, values of Nk{N) can be arbitrary numbers be- 
tween and 1; however, if we switch the interaction off 
{V=0), we revert to the single-particle model (^) in which 
iVfe=0, 0.5, or 1. 

From the shift-symmetry (^), it is obvious that the 
change in the ground-state energy, occurring when the 
interaction is switched on, is the same for all even sys- 
tems (Ceven) and the same for all odd systems (Codd) and 
hence 

S3p(7V)-t-Cevc„ for N^2n, 
Acq(iV) = { (42) 
Bsp{N) + Codd for N^2n+1, 



where Bsp(N) is given by Eq. (|2|). 
Filters (||), (^), and (||) now give 

A|2(iV) = Ag)(iV)-fCodd-Ce 



hence 



Aicq(iV) = Codd 

Set^iN) = <5eieq(iV) = d. 



(43) 



(44a) 
(44b) 



It is quite remarkable that the above argumentation does 
not at all depend on details of the two-body interaction, 
and that - irrespective of the interaction - filt ers ( 5|) and 
(^ correctly separate the i ntera ction effects (44a) from 
the single-particle spacings (44b). 

This generic result does not depend on whether any ap- 
proximations are used to obtain the ground-state energies 
of interacting systems. In particular, the BCS mean-field 
results (obtained for interaction V) will also obey the pat- 
tern presented in Eqs. (Q) exactly; only the value of the 
interaction-energy difference Codd — Ccvon may be differ- 
ent from the exact result. 

Similarly, the BCS results obtained for the seniority- 
pairing interaction ( |ll|) also follow the same pattern; 
however, in this case we have to additionally ensure that 
the phase space in which the pairing correlations are al- 
lowed to develop is the same for all particle numbers. 
Indeed, one usually solves the BCS equations in a finite 
phase-space window, adjusting the interaction strength 
C to that window. For the results (^) to be valid, we 
have to always use the same number of levels below and 
above the Fermi level, independently of the number of 
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particles. Actually, such a prescription is often used in 
realistic BCS or HFB nuclear structure calculations. 

As discussed above, in the case of the equidistant-level 
model, the energy-filter Se reproduces the value of the 
level spacing d regardless of the detailed structure of V 
and details of the many-body approximation used. The 
reason for this is the shift-symmetry of the Hamiltonian. 
Of course, if this symmetry is broken either by assum- 
ing the finite Hilbert space (see the following section) 
or by introducing the explicit symmetry- violating terms, 
one cannot a priori expect the above conclusion to hold 
exactly. 



It is also worth noting that the OES (44a) is defined in 
terms of the pairing (interaction) energy in the even and 
odd system, and not in terms of the HFB or BCS pairing 
gap A' (defined as half of the lowest excitation energy). 
Only for strong pairing correlations are the pairing gaps 
and the values of the OES the same (see Ref. ||ll|l). In 
nuclei, pairing correlations never reach such a limit, and 
the OES and BCS pairing gap in an odd system can be 
significantly different. 



D. Deformed shell (finite non-degenerate spectrum) 

In this section, we investigate the deformed-shell- 
plus-pairing Hamiltonian that contains the single-particle 
term (|lj) (for twofold degenerate single-particle states 
available for the pair scattering) and seniority-pairing in- 
teraction (|ll|). In this case, the analytic solution does not 
exist, but the exact eigenstates can be found numerically 
using, for example, the Richardson method [2S[| or by 
performing a direct diagonalization. The latter approach 
cannot be applied when the number of single-particle lev- 
els is large, because the dimension of the Hilbert space 

grows as ^^^^2^ ' ^'^^ ^^^^ puts the practical limit at 

r2 ^ 20. imaginary-time method 

A particular variant of this model is when the single- 
particle spectrum is uniformly spaced (finite-number 
equidistant-level model), i.e., = kd (fc=l,...,f2). In 
this work, we are mainly interested in the ability of 
the energy-spacing filter (^) to extract the single-particle 
spectrum from the total binding energies. Therefore, be- 
low, we also study the case of a nearly-equidistant spec- 
trum, in which Ck = kd except for the seventh level 
shifted up in energy by d/i (i.e, ey = 7.25d). The value 
of the single-particle spacing d constitutes a convenient 
energy scale, and below, the results will be expressed as 
ratios of all energies and parameters with respect to d. 

Figures ^ and g show the results of the exact calcula- 
tions for the nearly-equidistant spectrum of Q—16 levels. 
The behavior of the three-mass filter (|l|) is illustrated in 



Fig. |4|, where values of A^j^j are compared to 



which is the "equivalent" pairing gap. (In the BCS limit, 
A" becomes the gap parameter.) 

For G/d=0.1 (weak pairing ), values of A^;^f(2n -t- 1) 

nicely follow the low values of A", while those of A^g£(2n) 
are clearly influenced by the single-particle spectrum and 
do not at all reflect the smallness of the pairing correla- 
tions. In particular, the fluctuation around A^=14 {n—7) 
clearly shows up in A^^f(2n) and is absent in A^'^f(2n-|-1). 
In the case of intermediate pairing {G/d=0.3), 



^dct(2" + 1) behaves rather smoothly, while A|j^j(2n) 
zigzags in the region of irregularity in the spectrum. No 

direct correspondence between the values of A^^j and 
A*^ can be found here. However, with increasing val- 

(3) 

ues of G, i.e., when the static pairing sets in, Ajj,^ closely 
approaches A". This is nicely illustrated for G/d=0.5. 
Only in the case of relatively strong pairing correlations 
(but still far from the strong pairing limit discussed in 

I k f3) 

Sec. |III A 2D is the fluctuation in A^^j barely visible. 



In the case of weak pairing correlations, the corre- 
sponding energy-spacing filters (^), shown in Fig. H, very 
well reflect the structure of Hsp- With an increasing pair- 
ing strength, the pairing Hamiltonian gives rise to a dif- 
fused Fermi surface (i.e., it smears out the single-particle 
occupations). Consequently, the information about the 
details of the single-particle distribution is then expected 
to be washed out. This is clearly seen in Fig. ||, where the 
nearly-equidistant single-particle spacings are marked by 
dots. 

However, it is seen that even in the case of relatively 
strong pairing correlations, (G/(i=G.4), the symmetric 
energy-spacing indicator Se gives a qualitative descrip- 
tion of the spectrum. Even at G/d=0.5, the zigzag in 
Se appears in the right place. As far as the asymmetric 
filters Se^ are concerned, their behavior is more strongly 
influenced by G, and the resulting particle number de- 
pendence may make them less useful measures of the 
spectral properties. 

In order to assess the quality of the quasi-particle 
(BCS) approximation for the binding-energy indicators, 
in the same model we have also carried out the BCS cal- 



culations. As seen in Fig. 



6, it is only at large values 



of G that the A^^^ filters (1) which are applied to the 
BCS energies approach the exact results. Especially at 
intermediate values of pairing strength, where the static 
pairing vanishes in odd-iV nuclei, BCS becomes a rather 
poor approximation. However, it is clear that values of 
(2n) are affected by the single-particle spectrum, 



,(3) 

dof-l-BCS 



while those of A 
it. 



(3) 

dof-l-BCS 



(2n -|- 1) are rather insensitive to 



The agreement is significantly better for the energy- 
spacing filter Se (Fig. 0). Again, as discussed above, at 
large values of pairing strength, Se is only a qualitative 
measure of the single-particle splitting. 



A"^GV(PtP)-[7V/2], 



(45) 
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1. Second-order perturbation results 



Following Sec. [II A, one can derive analytic expres- 



sions for the binding energy filters which should be valid 
in the limiting cases of weak and strong pairing. Partic- 
ularly simple are the weak-pairing expressions. Namely, 
Eqs. (|l|) yield: 



A!,ek2n+1) 



^G+^[S{n) + S{n~n-l)], (46a) 



def 



(2n) 



d+G 
2 ^ 4d 



[S{n) + S{n - n)] , (46b) 



where 



n 



(47) 



Using Eq. (|T|), the single-particle splitting filters are 

G2 



Se+ = d + 
6e^d + 



2d{n - n) ' 
2dn 

G^ 

2dn{n - n) ' 



(48a) 

(48b) 
(48c) 



Note that in this case, the single-particle splitting filters 
( p8| ) are always greater than the original spacing d. In 
particular, Se~^ increases with n, Se~ decreases with n, 
and Se has a minimum in the middle of the shell. 

In the strong-pairing limit of the equidistant-level 

(3) 

model, the zero-order expressions for A^^f and 6e are 



given by Eqs. (^3|) and (^4|). The corresponding higher- 
order corrections can be easily derived using Eqs. (^l|)- 
( pO| ) after noting that 



n + i 



-d and err 



12 



(49) 



Figure [s] shows the accuracy of the second-order ex- 
pressions for the binding energies in the equidistant-level 
model with f2=16 in the weak-pairing (left) and strong- 
pairing (right) limits. It is seen that, for large numbers 
of particles, the low-rj expansion is fairly accurate even 
for ?7~0.4. Also the r]~^ expansion seems to work very 
well even for large values of d/G. different vertical 



E. Pairing-plus-quadrupole model 



In the previous sections, we have considered models 



having either de gener ate (Sec. IIIB) or arbitrarily fixed 
(Sees. [lie and [II D| ) single- particle spectra. However, 
in real systems, the single-particle energies do depend on 
the numbers of particles, and the energy-spacing filters 
(a) extract from masses the single-particle energies which 



include the particle-number dependence, i.e, at every A'', 
the distance to the next available single-particle level is 
obtained. 

In order to account for such effects, in this section we 
analyze the results of the exact diagonalization of the 
pairing-plus-quadrupole (PPQ) Hamiltonian, 



-GP^P - kQ ■ Q, 



(50) 



in a single-j shell. The PPQ model is well known and 
has been used many times to test collective properties 
of fermion systems; here we use the version defined in 
detail in Ref. Calculations were performed in the 

j=19/2 shell (ri^lO). The strength of the quadrupole- 
quadrupole (QQ) interaction k provides here a suitable 
energy scale, and in the following all the energies and 
parameters are expressed as ratios with respect to k. 

Exact ground-state energies for all particle num- 
bers N=0,l,...,20 and G=0 are plotted in Fig. |. 
All the even systems have the ground-state spins of 
/g.s=0, while in the odd systems, the ground-state spins 
are Ig,s^\N—n\/2+n/2; this corresponds to rotational 
bands based on oblate sequences of deformed single- 
particle levels (see the discussion in Ref. [p0[ ). 

Even in a large energy scale of Fig. Bl, the effect of 
the twofold Kramers degeneracy and the OES effect are 
clearly seen. Note that in this model the pure QQ in- 
teraction generates a (weak) OES effect. Indeed, in a 
relatively small phase space, the QQ interaction has a 
tangible pairing component. 

In order to quantify the OES and mean-field effects. 
Fig. 1^ shows values of the three-mass filter (^ obtained 
for different pairing strengths G. Ground-state energies 

(3) 

of even and odd systems are used to calculate AppQ for 
each value of G. For even systems, the ground-state spins 
equal for all values of A^ and G, while for odd systems, 
the values of /g.s=19/2 (for all A^) replace at large G 
the values of Ig,s=\N—n\/2+n/2, which characterize the 
G=0 solutions. This corresponds to a gradual transition 
from deformed to spherical shapes and from rotational 
to seniority-like excitation spectra. 

In Fig. |l^, one can clearly see that the OES increases 
almost linearly with G, while the pattern of alternating 

(3) 

larger and smaller values of Appq is almost independent 

of G. Both these features of Appg are explicated by using 
filters (H) and ( |6b|) which give values of Appq and Seppq 
plotted in Fig. |11|. One can very well see the almost 
linear dependence of the OES on the pairing strength G, 
Fig. Ul^(a), and a very weak G-dependence (apart from 
N—10) of the single-particle energy spacings, Fig. [ll|(b). 

The PPQ model exhibits several features pertaining 
to two kinds of phase transitions. First, the static pair- 
ing correlations set in at critical values of the pairing 
strength G. Depending on the number of particles, this 
phase transition occurs at about G/k=0.03-0.05 for even 
particle number. Second, the transition from deformed 
to spherical shapes occurs at slightly higher values of G, 
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i.e., at about G/k=0.08, 0.10, 0.12, and 0.12 for N=4:, 
6, 8, and 10, respectively. In addition, deformed equilib- 
rium shapes are oblate for N—2, 4, 6, and 8, and prolate 
for 7V=10. (Note that apart from a linear dependence 
of energies on N, the PPQ Hamiltonian is exactly sym- 
metric with respect to the particle-hole transformation 
N^2n-N |l|.) 

Since the system is finite (and fairly small for that 
matter), the phase transitions are hardly visible in the 
exact results of Fig. |ll|. However, when filters (^ and 
( |6b| ) are applied to the mean-field (HFB) ground-state 
energies (Fig. 12), the phase transitions become visible 
as sudden increases in Appq+hfb (pairing transition), 
and degeneracies of ^eppq+HFB (shape transition). 

As one can see, a comparison of the exact and HFB 
ground-state energies (Figs. |l^ and 12) is very instruc- 
tive when it is based on comparing the corresponding 
filters (||) and (6b). It turns out that in the PPQ model 
the HFB method reproduces quite well the OES and 
the single-particle properties simultaneously. Some de- 
viations occur only near the phase transitions, where it 
is well known that the mean-field approximation is not 
accurate. 

When analyzing exact solutions for systems interacting 
with two-body interactions, or when similarly analyzing 
the experimental data, one does not have a priori ac- 
cess to the single-particle energies or to the single-particle 
energy-spacings. In fact, the single-particle energies are 
concepts that appear naturally only in the mean-field ap- 
proximation. Therefore, in order to assess the meaning 
of numbers obtained from the energy-spacing filter (|6b|), 
one should compare them with the energy-spacings cal- 
culated directly from the mean-field spectra, i.e., with 
the differences 

SeiN) - £jv/2+i(iV) - eN/2{N), (51) 

[cf. Eq. (0)] where ek{N) is the fcth (twofold degenerate) 
single-particle energy obtained within the mean-field ap- 
proximation. Since a consistent application of the mean- 
field approximation to an odd system requires the time- 
reversal breaking, the Kramers degeneracy is lifted in odd 
systems. Consequently, in Eq. ( pl[ ) one should only use 
the single-particle energies obtained self-consistently for 
systems with an even number of particles. 

In Fig. |l^(b) we show the differences ( |5l| ) calculated 
from the single-particle spectra of canonical HFB ener- 
gies obtained in the PPQ model. (In fact, since all of 
the PPQ-fHFB equilibrium solutions conserve the axial 
symmetry, the PPQ+HFB method reduces to the sim- 
ple BCS approximation, and the canonical energies are 
equal to the eigenenergies of the mean- field Hamiltonian.) 
It is clear that the energy-spacing filters (applied either 
to the exact or to the HFB total energies) give results 
similar to the differences of canonical energies only for 
deformed shapes. Whenever the mean-field solutions be- 
come spherical, differences ( ^l| ) collapse to zero, as ex- 
pected, while the energy-spacing filters give non-zero val- 
ues. This result is easy to understand: in the spherical 



limit, the results of the PPQ model should resemble those 
of the seniority model. While the single-particle energies 
are degenerate (the Hilbert space consists of one j-shell 
only), the energy-spacing filter sho uld b e proportional to 
the pairing strength G (see Sec. IIIB). On the other 
hand, before the transition to sphericity, the results pre- 
sented in Figs. |l](b), |2|(b), and |3|(b) are encouragingly 
similar. Note that deviations obtained at iV=10 should 
again be attributed to the differences in predicted shapes 
(prolate vs. oblate). 

Finally, in Fig. |l^(a) we show the values of the 
HFB "equivalent" gap parameters (|5|) calculated in the 
PPQ-I-HFB model. It is seen that A° significantly un- 
derestimates the magnitude of the OES effect, and more- 
over, it exhibits some particle- number dependence which 
is absent in the exact results. 

In light of the above discussion, the very weak aver- 
age dependence of Appq on N (except for the OES, 
of course), shown in Fig. |l^, can be given by a very 



(3) 
PPQ 



simple explanation. In the weak pairing limit, A 
is small, and its overall particle-number dependence is 
much weaker than the even-odd effect (see Fig. P). On 

(3) 

the other hand, in the limit of strong pairing, Aj3pQ is 

expected to approach the seniority limit in which A^^-* 
depends only on the number parity but not on N. Note 
that the results of the non-degenerate model shown in 
Fig. ^are very far from the spherical limit since G/d < 1. 



IV. CONCLUSIONS 

This work contains the analysis of an interplay be- 
tween pairing and mean-field effects on binding energies 
of many-fermion systems. While most of our discussion 
is concerned with nuclear systems, the main conclusions 
also apply to other finite-size superconductors such as 
grains. 

The analysis of binding energies of several exactly- 
solvable Hamiltonians (allowing variations in the magni- 
tude of pairing correlations) demonstrates that the three- 
mass indicator, A'^^^(27i + 1), is indeed an excellent mea- 
sure of pairing correlations, and the symmetric filter Se 
adequately extracts the effective single-particle spacings 
from the measured binding energies of deformed nuclei. 
According to the analysis of nuclear masses, the mean- 
field contribution to the OES is significant in light nuclei, 
but it is reduced to ~100-200keV in heavy deformed 
nuclei due to the relatively close spacing of single- 
particle levels. 

For deformed nuclei with weak and intermediate pair- 
ing correlations, there is a nice consistency between the 
approximate results obtained within the mean-field ap- 
proach (BCS or HFB) and the exact results. However, for 
strongly paired nearly-spherical nuclei, there is no clear 
correlation between A('^)(2n -I- 1) and the pairing defor- 
mation A°. Also, the energy-spacing filters are superior 
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over the mean-field single-particle energies in assessing 
the single-particle properties of the system. 

Approximate expressions of binding-energy indicators 
have been derived in the limits of weak and strong pair- 
ing. These formulae nicely explain the gross particle 
number dependence seen in the exact results. 
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FIG. 1. Neutron OES filter A^^p (|) plotted as a function 
of neutron number N (A-odd). Results for different isotones 
are marked by dots. The average values of Aoxp are indicated 
by gray bars. Experimental data were taken from Ref. ||2^. 
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FIG. 2. Same as in Fig. ^ except for the proton OES filter 
^exp plotted as a function of Z [Z-odd). 
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FIG. 4. Three-mass fihers A^^^', ^ (thick solid and dot- 
ted lines) calculated for the exact binding energies in the 
deformed-shell-plus-pairing model with 0=16 for the case of 
weak (G/d=0.1), intermediate {G/d=0.3), and strong pair- 
ing (G/d=0.5). The single-particle spectrum is uniform 
(efe — dk), except for the seventh level which is shifted up 
in energy by d/4 (i.e, e? = 7.25d). The equivalent gap pa- 
rameters (lisn are shown by thin lines. 
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FIG. 3. Average values of neutron and proton OES filters 
(^) plotted as functions of the neutron and proton numbers 
(both and Z odd), respectively. Experimental masses were 
taken from Ref. ifesl. 



FIG. 5. Similar to Fig. ^ except for the energy-spacing 
fihers (|) 5edef calculated for G/d=0.1, 0.2, 0.3, and 0.4. 
The nearly-equidistant single-particle spacings are marked by 
dots. 
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FIG. 6. Three-mass filters A^^j. (|) calculated for the 
exact (thick lines) and BCS (thin lines) binding energies 
in the deformed-shell-plus-pairing model with f2=16 and a 
nearly-equidistant spectrum. Solid and dotted lines show re- 
sults for even and odd values of A*', respectively. 



FIG. 8. Relative error (in percent) of the second-order 
expressions for the binding energies of the equi distant level 
model with S7=16 in the w eak (lef t, G/d<^l; Sec. [II A 1) and 
strong (right, d/G<^l; Sec. [II A 2) pairing limits. In the weak 
pairing limit, calculations were performed for A'^=4, 5, 14, and 
15. To realize the strong-pairing situation, only large particle 
numbers, A'^=14 and 15, were considered in the d/G<^l case. 
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FIG. 7. Similar to Fig. | except for the energy-spacing fil- 
ters (^ ^edef. The single-particle spacings are marked by 
dots. 
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FIG. 9. Exact binding energies Bppq (solid line) of par- 
ticles in the j=19/2 single-j shell interacting with the pure 
QQ interaction (G—Q). Energies obtained within the HFB 
approximation are shown with the dashed line. 
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FIG. 10. Exact values of the three-mass fiher A 
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the j = 19/2 PPQ calculated for the pairing strengths G/k, 
indicated at the right-hand side. 
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FIG. 12. Same as in Fig. |l^ except for the HFB results. 



FIG. 11. Exact values of the OES Appq (|) (a) and of the 
energy spacing 5eppQ ( |6b|) (b) calculated within the j=19/2 
PPQ model. 
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FIG. 13. Same as in Fig. |Tl] except for the Appq+hfb HFB 
order parameter ( ^ ) and the differences of HFB canonical 
energies (5£ppq+hfb 
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